Four Years of Extreme Ultraviolet Observations of Markarian 421. 

I: Spectral Analysis 
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ABSTRACT 

We analyzed the ~ 950 ks of spectroscopic data accumulated by the Extreme Ultraviolet 
Explorer (EUVE) satellite between 1994 and 1997 for the BL Lacertae object Markarian 421. 
The EUV spectrum is well detected in the 70-110 A (112-177 eV) range and can be fitted by 
a power law model plus an absorption feature in the ~ 71 — 75 A range. Previous studies of 
EUV absorption features in Mrk 421 and in the other EUV bright BL Lac object, PKS2155-304, 
explain this absorption feature as a superposition of Doppler-smeared absorption lines (mainly 
L- and M-shell transitions of Mg and Ne) originating in high-velocity gas clouds ionized by the 
beamed continuum of the associated relativistic jet. We show that, for example, Fe IX L could 
also be a possibility consistent with the marginal detection of oxygen absorption lines in the X- 
ray range. However physical models are highly sensitive to the assumptions on the photoionizing 
continuum and the surrounding gas. 

Subject headings: galaxies: nuclei - galaxies: active - galaxies: individual (Mrk421) 



1. Introduction 

Evidence for absorption features in the X-ray 
spectra of BL Lac objects at ~ 0.6 keV 1 has ex- 
isted for many years. This is a tantalizing situ- 
ation: BL Lac objects were originally defined by 
their lack of the sharp spectroscopic features that 
provide the normal tools for understanding astro- 
phyisical objects. Without this tool, only large 
multisatellite, multitelescope campaigns spanning 
20 decades of frequency have been able to con- 
strain the BL Lac emission mechanism. Absorp- 
tion features can provide a unique tool to study 
the presence and the dynamical properties of any 
gas surrounding the nucleus of the otherwise 'fea- 
tureless' BL Lac objects. The extremely large 
Extreme Ultraviolet Explorer {EUVE) satellite 



J PKS 2155-304 (Canizares & Kruper 1984), H1426+428 
(Sambruna et al. 1997) and PKS 0548-322 (Sambruna & 
Mushotzky 1998) 



database on Mrk 421 offers a so far a unique op- 
portunity to search for and characterize any sharp 
spectroscopic feature. 

Mrk 421 is one of the brightest, best known 
and closest BL Lac objects (z=0.0308, Ulrich et 
al. 1975). Mrk 421 was in fact the first BL Lac 
object found to have X-ray emission (Richetts et 
al. 1976) and its low redshift makes this object 
one of the few extragalactic sources detected at 7- 
rays above 300 GeV up to TeV energies (Punch et 
al. 1992). Mrk 421 has been extensively observed 
at radio (e.g. Zhang & Baath 1990), UV/optical 
(Maza, Martin & Angel 1978; Mufson et al., 1990) 
and X-ray frequencies (e.g. Mushotzky et al., 
1979; George, Warwick & Bromage 1988, Taka- 
hashi et al. 1996, Guainazzi et al. 1999). Mrk 421 
shows optical polarization, a flat radio spectrum 
and significant time variability, characteristics of 
the blazar class (Makino et al., 1987) and its par- 
ent galaxy has been identified as a giant ellipti- 
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cal (Ulrich et al. 1975; Mufson, Hutter & Kondo 
1989). 

In the EUV band (~ 60 - 750 A) Mrk 421 was 
detected during the ROSATWFG survey (Pounds 
ct al. 1993, Pye et al. 1995) and the EUVE all- 
sky survey (Marshall, Fruscione & Carone 1995, 
Fruscione 1996). It is one of the strongest EUV 
extragalactic sources and was immediately recog- 
nized as a prime target for subsequent EUV spec- 
troscopic studies. For 4 almost consecutive years 
(1994, 1995, 1996 and 1998) Mrk 421 was the tar- 
get of multiwavelength campaigns covering radio 
to TeV energies (Macomb et al., 1995, Kerrick et 
al. 1995, Buckley et al. 1996, Takahashi ct al. 
1999) and the EUVE satellite (which gives simul- 
taneous photometric and spectroscopic data) was 
always part of the campaigns. In addition, in 1995 
and in 1997 Mrk 421 was observed serendipitously 
with the EUVE photometers (the "scanners"). 
The large amount of data accumulated for Mrk 421 
by the EUVE satellite to date (~ 1100 ks of public 
imaging data and ~ 950 ks of public spectroscopic 
data) represents the best coverage for any BL Lac 
object at high energies. Less than half of this data 
has been analysed to date (Fruscione et al. 1996, 
Kartje et al. 1997, hereafter K97) revealing that 
(i) the source exhibits a flare-like behavior (ii) the 
strong EUV variability is correlated with soft X- 
ray and TeV-energy emission, and (iii) the spec- 
trum may have strong absorption features between 
65 A and 75 A (0.16-0.19 keV) (K97). However in 
the EUV, all the previous analysis of BL Lacs can 
only infer the presence of absorption from the ex- 
trapolation of the X-ray spectrum (e.g. Konigl et 
al 1995 for PKS2155-304 and K97 for MRK 421). 

The purpose of this paper is to analyze, in an 
homogeneous way, all of the 950 ks of EUVE pub- 
lic spectroscopic data for Mrk 421 . The results 
of this paper will also provide a useful reference 
for the calibration of other satellites observing in 
this energy range. In particular Mrk 421 is a cal- 
ibration target for both the Chandra and XMM 
satellites, whose Low Energy Transmission Grat- 
ing (LETG) and Reflection Grating Spectrometer 
(RGS) partially overlap with the EUVE energy 
band. 

In Section 2 we present all the EUVE spectra of 
Mrk 421 taken from 1994 to 1997; in Section 3 we 
discuss the results and in Section 4 we give a brief 
summary and our conclusions. All the photomet- 



ric data, the corresponding lightcurves and a de- 
tailed variability analysis are presented in a com- 
panion paper (Cagnoni, Papadakis & Fruscione, 
2000, hereafter paper n). 

2. Observations 

The EUVE satellite was launched in June 1992 
and is still acquiring astronomical data in the 
wavelength range from ~ 60 to 750 A. Onboard 
EUVE there are 4 telescopes: three of them are 
co-aligned photometers observing in 4 different 
EUV bandpasses, while the fourth is the Deep 
Survey/Spectrometer (DS/S, see e.g. Welsh et al. 
1990) mounted orthogonally to the scanners. (Fig- 
ure 1 of paper II shows the DS effective area com- 
pared to the scanner effective area) . In addition to 
a photometer (10 % bandpass 68-178 A) the DS/S 
telescope is also equipped with three spectrome- 
ters (Hettrick & Bowyer, 1983; Abbott et al. 1996) 
covering the "short" (SW: 70-190 A), "medium" 
(MW: 140 - 380 A) and "long" (LW: 280 - 760A) 
EUV wavelengths. This configuration allows si- 
multaneous imaging and spectroscopy with a spa- 
tial resolution of <~ 1' and a spectral resolution of 
A/AA ~ 200 at the short wavelengths, i.e. 0.35 A 
at 70 A. 

Further details on the science instrumentation and 
performances can be found in Malina (1994), Sirk 
et al. (1997) and Abbott et al. (1996). 

The scanners were used to carry out the EUVE 
all-sky survey (July 21 1992 - January 23 1993) 
and are currently utilized in the Right Angle Pro- 
gram (RAP), that is the observation of targets of 
scientific interest with the scanner telescopes while 
the DS/S telescope is conducting the primary sci- 
ence observation. The DS/S has been used for 
a deep EUV survey along the ecliptic during the 
all-sky survey and for pointed spectroscopic and 
imaging observations since then. 

We use throughout this paper all of the pub- 
licly available DS/S data on MRK 421: from 1994 
to 1997, a total of 950 ks. We are not consid- 
ering the most recent observation (April 1998) 
since the data is presently being analysed by the 
original proposers (Marshall et al. in preparation) . 
Mrk 421 was observed by EUVE several times 
from 1994 to 1997, four times with the DS/S and 
twice with the scanners; Table 1 summarizes the 
observations. 
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Because of strong absorption by the interstellar 
medium (ISM) along the line of sight (even at 
the low N H = 1-45 x 10 20 cm" 2 Elvis, Wilkes 
& Lockman 1989), Mrk 421 was detected only 
at the shortest wavelengths: in the Lex/B filter 
of the scanners (10% bandpass 58-174 A), in the 
Lex/B filter of the DS instrument (10% bandpass 
67-178 A) and, relevant for this paper, in the SW 
spectrometer (70-190 A) 2 . 

2.1. Energy spectra 

After rejecting the time intervals with high par- 
ticle background, such as the satellite passages 
over the South Atlantic Anomaly, and correcting 
the data for instrumental dcadtime and telemetry 
saturation, we measured a total effective exposure 
of 315858 s, 506780 s, and 354806 s for the 1994, 
1995 and 1996 observations respectively (note that 
the first 1995 and the 1997 observations were per- 
formed with the scanners, instrument with imag- 
ing capabilities only). 

We extracted all the spectra in a homogeneous 
way from the two-dimensional SW detector image, 
using the EUVE Guest Observer Center software 
(IRAF/EUV package) and other standard spectro- 
scopic IRAF tasks. The spectrum was extracted in 
a 14 pixel wide aperture and the background ob- 
tained by linearly fitting an average background 
in two regions (~ 80 pixel wide) of the detector, 
one on each side of the spectrum. The resulting 
wavelength-calibrated spectra (Figures 1, 4 and 8) 
were convolved with the proper effective area and 
binned over 1 A. For the 1994 on-axis observation 
(Figure 1) we used the on-axis effective area de- 
rived from in-orbit calibrations for the 75 A- 110 A 
range. The 1995 and 1996 observations were per- 
formed off-axis (a technique that slightly increases 
the short- wavelength coverage, down to ~ 70 A). 
This results in a different effective area (Marshall 
et al. 1999) that we applied here for the first time. 
The lcr error bars on the flux measurements shown 
in the figures of all the spectra were computed us- 



ing the formula 



2 A graphical representation of the simulated DS/S and scan- 
ner responses is presented in Figure 2 of paper II, where we 
convolved the Mrk 421 energy spectrum corrected for the 
ISM absorption with the DS/S and scanners effective area. 
No signal is present at wavelengths longer than 



110 A 



a(X) 



^[S x + B x /(1 + H sp /H B )} 

AeffW 



where S is the signal from the target, B the aver- 
aged background, both in counts/s, A e ff the SW 
effective area in cm 2 and H sp and Hb are respec- 
tively the width in pixels of the extraction region 
for spectrum (14 pixels) and background (160 pix- 
els) in the direction perpendicular to the disper- 
sion. 

We modeled the EUV spectrum of Mrk 421 
using version 2.0 of the SHERPA modelling and 
fitting software (recently developed at the Chan- 
dra X-ray Center). We initially used an absorbed 
power law of the form 



f(A) = f(A )A«- 1 cxp{-[I]xN(X) ( Tx)} 



photons cm 2 s 2 A 



where f(Ao) is a normalization factor. We used 
A = 80 A (E =0.155 keV). a is the energy index 
and the absorption is characterized by a column 
density N(X) and an absorption cross section ax 
for each element. We included in the model H, He 
I and He II (Rumph, Bowyer and Vennes, 1994) 
and heavier elements (Morrison and McCammon, 
1983). We fixed the Galactic hydrogen column 
density at N H = 1.45 x 10 20 cm" 2 (Elvis, Wilkes 
& Lockman 1989) and the ratios Njj e i/ Njfi=0.1 
and Nfj eii / Nui=Q. 01. The best fit values for the 
simple absorbed power law model are summarized 
in Table 2 together with their la uncertainties for 
1 interesting parameter. Two-parameter x 2 con- 
tour confidence levels at Xmm+1' +4.61 and +5.99 
- corresponding respectively to the la uncertainty 
for each single parameter and to the 90% and 95% 
joint confidence levels - are shown in Figure 2, 5 
and 9 for the 1994, 1995 and 1996 data sets re- 
spectively. 

2.1.1. 1994 energy spectrum 

Our best fit flux in 1994 is slightly higher than 
the previous analysis of Fruscione et al., 1996 

( f soA = 6 ' 44 ± °' 33 X 10 ~ 3 P hotons cm ~ 2 s_1 
A -1 compared to their f ^ = 290 ±15 /iJy= 

5.47±0.28 x 10~ 3 photons cur 2 s" 1 A" 1 ) but our 
slope is within the errors (a = 0.85+Q 72 compared 
to their a = 1.4 ±0.8). The la error bars reported 
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in Fruscione et al., 1996 were not corrected for the 
effective area; the small difference in the value of 
the normalization can be possibly explained by the 
difference in the errors and a possible difference in 
the wavelength range used to model the spectrum. 
This is the only spectrum for which the addition 
of an absorbed interved gaussian line does not im- 
prove the fit (Tab. 3). 

2.1.2. 1995 energy spectrum 

The 1995 EUVE DS/S observation was part of 
a large multiwavelength campaign during which a 
powerful TeV flare occured (Buckley et al. 1996). 
Simultaneously an X-ray flare developed (Taka- 
hashi et al. 1996), and the EUV followed with a 
short delay (Buckley et al. 1996). As illustrated 
in figure 3, the EUVE observation spanned the 
flare (~ 104 ks), the decay (~ 218 ks) and the 
following quiescent state (~ 185 ks) The obser- 
vation was performed ~ 0.3° off-axis and as ex- 
pected the signal reaches down ~ 7fjA. This spec- 
tral region is crucial since it covers the location of 
the wide EUV absorption feature proposed by K97 
using less than half (the flare decay) of the 1995 
data. In our analysis we corrected the data with 
the ~ 0.3° off-axis effective area (Marshall et al. 
1999) and the resulting spectrum, binned over 1 
A, is presented in Figure 4. 

In order to compare our data analysis to that of 
K97, and in an attempt to look for changes dur- 
ing the flare, we divided the 1995 observation into 
three parts (Figure 3): 

Flare: April 25-28 (~ 104ks), the time interval of 
the EUVE flare; 

Decay: April 29 - May 6 (~ 218 ks), the time in- 
terval analysed by K97; 

Quiescence: May 7-13 (~ 185 ks), the remaining 
data. The EUVE spectra for the 3 parts are shown 
in Figiure 6 a,b and c. 

Table 1 lists the exposures times and Table 2 
summarizes the best-fit parameters for each one 
of the three parts. In Figure 7 we plot the two- 
parameter x 2 contour confidence levels. 

2.1.3. Comparison with previous results 

The 'decay' section of the 1995 spectrum de- 
serves comparison with the analysis of K97, given 
the importance of the absorption line evident in 
their Figure 2. We show the effect of using the 



improper ~ 0.3° off-axis effective area in Figure 
6b; we compared the EUVE spectrum obtained 
using the on-axis effective area (open dries) as 
in K97, and the 0.3° off-axis effective area (filled 
squares). The effect of using the appropriate cal- 
ibration is to increase the flux below ~ 80 A and 
so weaken the absorption feature (~ 45%). An- 
other effect to consider is that K97 evaluated the 
level of the continuum in the EUV range by adopt- 
ing an «euvx (~ 1-32) slope derived from the 0.15 
kev EUV flux and the 1.5 keV X-ray (ASCA) flux, 
which resulted in a strong absorption feature (see 
also §3 for further discussion on this point). 

Our best fit slope a = 2.19 ± 0.5 (Table 2) is 
flatter than the a = 3.5±0.8 in K97. If we consider 
the same wavelength range (70-95 A instead of 70- 
110 A) we find a = 2.54±0.56, consistent with K97 
result. 

A simultaneous fit of an absorbed power law 
plus an absorbed (inverted) gaussian should pro- 
vide a more reliable result. If the position of 
the gaussian is a free parameter then the reduced 
X 2 goes to 1.17 and a narrow (FWHM=1.8 ± 0.5 
A) absorption feature is found at ~ 72.7 ± 0.3 
A while the slope becomes significanlty flatter 
(a = 1.62 ± 0.59) (see Tab.3). This is the first 
time that an absorption feature in the EUV range 
is derived from EUV data alone. 

Both for PKS 2155-304 (Konigl et al. 1995) and 
for Mrk 421 (K97) the presence of absorption fea- 
tures was previously derived on the basis of an 
extrapolation of the X-ray power law and not di- 
rectly from EUVE data. However there is strong 
evidence that the X-ray spectrum of the high- 
energy peaked BL Lac objects, including Mrk 421 
, is concave (Sambruna et al. 1994, 1997; Tashiro 
et al. 1995; Takahashi et al. 1996a; Giommi et al. 
1998; Guainazzi et al. 1999) and a gradual steep- 
ening with energy agrees with the "synchrotron 
self Compton" scenario, which has often been used 
to explain the spectral energy distribution of BL 
Lac objects (Ghiscllini, Maraschi & Treves 1985, 
Ghiscllini 1989). Therefore, a simple power law 
model is an unreliable extrapolation from X-rays 
to the EUV. 

In order to verify this assumption for Mrk 421 
, we reanalyze, as an example, one of the many 
one-hour-long 1995 ASCA observations (the one 
on 25-26 April) undertaken simultaneously with 
the EUVE observation (1995 'flare'). This in or- 
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der to investigate the spectral slope below 1.5 
keV (Takahashi et al. 1996b and 1999 only pub- 
lished data for 1.5 keV and above). First we fit- 
ted the ASCA spectra with an absorbed broken 
power law and the best fit requires a break posi- 
tion at Ebreak = 1-77 ± 0.07 keV and the energy 
slopes varies from 1.44 ± 0.16 above the break to 
0.98 ± 0.02 below the break (consistent with what 
found by Takahashi ct al.). The broken power law 
model however does not give a good fit of the spec- 
trum below 1.5 keV; a better result is found if we 
assume a "curved spectrum" parametrized as 

F(E) = E- {J{E)ai ° w+{1 - f{E))ah ^ h) 

where f(E) = (1 - exp(-E / Eo)) 13 , a low and a high 
arc the low and high energy asymptotic slopes 
and (5 is the curvature radius in the energy space. 
The best fit values for the curved model with 
fixed (3 — 1.0 (e.g. Guainazzi et al. 1999) are 
a low = 0.57 ± 0.06, a high = 1.34 ± 0.13 and E = 
1.78 ± 0.18 (see fig. 10); these values show that 
the slope flattens toward lower energies and that 
an extrapolation with a simple power law is not an 
accurate way to predict the EUV continuum level; 
the low energy asymptotic slope (ai ow — 0.57) is 
generally in agreement with the flat slope derived 
from the EUV spectrum (a = 0.71 ± 0.61 from 
Tab.3). 

Note that in the case of PKS 2155-304 there is 
an additional problem: the ROSAT spectral slope 
used for the extrapolation was not measured si- 
multcneously to the EUV spectrum and the spec- 
tral slope has been observed to harden when the 
source intensity increases (Giommi et al. 1998). 

2.1. 4- 1996 energy spectrum 

In 1996 Mrk 421 was again observed ~ 0.3° off- 
axis with a total effective exposure of 354806 s. 
The 1 A binned spectrum and the best fit pa- 
rameters model are plotted in Figure 8; the cor- 
responding contour confidence levels are shown in 
Figure 9. 

2.1.5. The EUV absorption feature 

We performed several fits using an absorbed 
power law model plus gaussian for all the observa- 
tions. The results are listed in Table 3. In every 
case, except for the 1994 observation, the addition 
of a gaussian absorption line improves the fit. The 



power law slope is in the range a = [—0.15, 1.62] 
which is not too distant from the slope measured 
from X-ray data, (see e.g. Guainazzi et al.with 
SAX data which gives a = 0.3 — 1.1 depending 
on the spectral model) . The best-fit center of the 
gaussian absorption lines is in the 71-75 Arangc 
with full width half maximum varying from 0.2 
to 6 A. We note however that the 6 A FWHM 
is relative to the global 1995 observation and it 
is probably due to the superposition of narrower 
slightly variable lines. 

3. Discussion 

3.1. The absorption feature 

Ever since Einstein Pioneering observations 
suggested the existence of an ubiquitous absorp- 
tion feature around ~ 0.55 keV in the X-ray spec- 
tra of BL Lacertae objects (Canizares & Krupcr 
1984, Madejski et al. 1991), much X-ray and EUV 
data has been collected with the goal of confirming 
the existence of such absorption features. 

3.1.1. EUV '/X-ray absorption in other BL Lac 
objects 

The first object to show evidence of an X-ray 
absorption feature was PKS 2155-304 (z = 0.117) 
in a 1980 Einstein grating spectrum (Canizares 
& Kruper 1984). The feature was detected at 
- 0.6 keV (A = 21 A) and extended for 50-100 
eV (A = 120 - 250 A); it was interpreted by the 
authors as O VIII Lya emitted by high veloc- 
ity material. Looking at a 1979 Einstein solid 
state spectrometer observation of the same object 
Madejski et al. 1991 found evidence for absorp- 
tion at ~ 0.55 keV, consistent with the grating 
result. 

Krolik et al. 1985 modelled the absorption as due 
to an outflowing subrelativistic wind beamed in 
our direction. 

PKS 2155-304 was subsequently observed by 
BBXRT in 1990 and the absorption feature at 
~ 0.55 keV was confirmed in those data (Made- 
jski et al. 1994), although instrumental effects 
cast some doubt on this result (Sambruna et al. 
1997, Weaver et al. 1995). 

In a long exposure (> 100 ks) Beppo-SAX spec- 
trum of PKS 2155-304 in the soft X-rays (Giommi 
ct al. 1998) the best fit model is improved (at 
99% confidence) when a notch at 0.55 ± 0.04 keV 
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is added to the model, although the magnitude of 
the feature is close to the instrument calibration 
limits. 

PKS 2155-304 was also observed by EUVE. A 
1992 30 ks observation gives hints of an absorp- 
tion feature at ~ 80 A (Fruscione et al. 1994). 
Two longer observations in June (130 ks) and in 
July (150 ks) 1993 suggested the existence of 2 to 
5 absorption features at a > 3a level in the range 
75 to 110 A (Konigl et al. 1995). Extrapolating 
the ROSAT X-ray spectral slope into the EUV 
the authors found evidence of a wide absorption 
feature from ~ 75 to ~ 85 A and interpreted it as 
the superposition of Doppler-smeared absorption 
lines (L-shell transitions of Mg and Ne as well 
as M-shell transitions of Fe) that originate in a 
high-velocity clumped outflow from the nucleus. 
According to this model the X-ray line spectrum 
should be dominated by an O VII Ka broadened 
and blueshifted feature. This corresponds to the 
feature observed in by Einstein and in 1990 by 
BBXRT, but the predictions are an order of mag- 
nitude greater than the observations; to match 
BBXRT measurement the oxygen must be present 
in subcosmic abundance (Konigl et al. 1995). 

Another BL Lac object to show an X-ray ab- 
sorption feature is H1426+428 (z = 0.129). Sam- 
bruna et al. (1997) reanalyzed a BBXRT 1990 
observation and found a broad feature around 
0.5 - 0.6 keV confirmed by an ASCA 1994 spec- 
trum, but not by a ROSAT 1993 spectrum. 

A third BL Lac, PKS 0548-322 (z = 0.069) 
(Sambruna & Mushotzky, 1998) confirmed the 
early discovery of an absorption feature from the 
Einstein spectrum (Madejski et al. 1991); this fea- 
ture is present in ASCA data and can be modeled 
with either an edge at ~ 0.66 keV and optical 
depth t <~ 0.3 or with a notch at ~ 0.82 keV and 
fixed width of AE = 0.1 keV and covering frac- 
tion f c ~ 0.1. The ASCA absorpton feature is 
at energies significantly higher than in previous 
X-ray observations, when the continuum was in 
a slightly higher state (by a factor ~ 1.4); this 
provides tentative evidence for variability in the 
absorption feature. 

3.1.2. EUV/X-ray absorption in Markarian J^21 

K97 interpret the absorption feature in the 
context of the Konigl et al. (1995) model for 
PKS 2155-304 as a superposition of Mg VIII and 



Mg IX absorption lines. The implied ionization 
state of the gas would require stronger O VII ab- 
sorption in the X-ray range than is allowed by 
observations and the authors invoke subcosmic 
oxygen abundance to explain its absence. 

Our analysis of the 1994, 1995 and 1996 EUVE 
spectra of Mrk 421 shows an absorption feature 
inferred directly from the EUVE data, (i.e. by 
modelling the EUV continuum and the absorption 
line simultaneously). The simplest model for this 
feature is a reverted Gaussian at about 71-75 A a 
maximum FWHM of the order of few A. 

In order to investigate likely absorption features 
in a BL Lac environment, we ran photoionization 
equilibrium models (Nicastro et al. 1999, Nicas- 
tro et al. 2000) which include the strongest 300 
emission and resonant absorption lines down to 
50 eV (oscillator strength > 0.1). We find that 
the predicted wavelengths and strengths of the 
absorption lines, both absolute and relative, are 
highly sensitive to the properties of the photoion- 
izing continuum and the absorbing gas (e.g. the 
ionization state and the assumed outflowing ve- 
locity). The number of transitions in the model 
also has a strong effect. However assuming the 
observed EUV-to-X-ray continuum and adjusting 
the gas paramters to maximize the relative abun- 
dance of OVII (~ 65%) both an OVII Ka absorp- 
tion line (EW=0.25A) and a FcIX L (A rcst = 82.43 
A, EW=0.2A) absorption line are expected respec- 
tively around 0.6 keV and 73 A (assuming an out- 
flow velocity of 0.15c). The equivalent width of 
the lines is too small to give significant detections, 
however, especially in the EUV regime, the pres- 
ence of many fainter satellite lines (not included 
in the model) could explain the observed feature 
For example between the Si VI line (A = 72.34 A) 
and the Mg VI complex (A 71.3 A), both included 
in the model, at least 15 additional lines from ion- 
ized Fe, Al, Be and Ne are not included (oscillator 
strength 0.01-0.1) and could be important espe- 
cially if the elements are abundant. 

Even though no strong absorption line seem to 
be visible in the X-ray spectrum, we note that 
both in the 1979 Einstein data (Madejski et al. 
1991) and in the 1997 and 1998 Beppo-SAX spec- 
trum (Guainazzi et al. 1999, Fossati et al. 2000) 
the fit improves when adding an absorption fea- 
ture centered at ~ 0.55 keV. 

In order to investigate any possible correlation 
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between spectral changes and source flux, we plot- 
ted in Fig. 11 the energy index as a function of to- 
tal EUV flux for the the power law plus gaussian 
models. However the large error bars prevent any 
statististically significant conclusion. 

4. Conclusions 

From EUVE data, we have found a clear ab- 
sorption feature at ~ 71 — 75 A in the spectrum of 
MRK 421. The absorption features detected in the 
EUV/X-ray spectra of BL Lac objects are often 
poorly determined statistically, as in the Beppo- 
SAX observations of PKS 2155-304 and MRK 421 
(Giommi et al. 1998 and Guainazzi et al. 1999), 
but it is interesting to notice that different instru- 
ments on different satellites point to the same pos- 
sible 0.55 keV feature and to a (probably linked) 
absorption feature in the EUV spectra of the same 
objects. 

We do not know if X-rays and EUV absorption 
features are present in all or many BL Lac objects, 
as initially suggested by Madejski et al. 1991, but 
to-date absorption features at 0.5 — 0.6 keV have 
been found - or at least statistically improves the 
spectral fitting - in four of the brightest BL Lacs 
objects (PKS 2155-304, H1426+428, PKS 0548- 
322 and MRK 421) and EUV absorption features 
has been found in two of them (PKS 2155-304 and 
MRK 421). Future observations covering the EUV 
to X-ray energy bands with much higher resolution 
and collecting with Chandra and XMM, 

will soon be able to answer this question and give 
more detailed observational data on which to de- 
velop new models. 
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Fig. 1— EUV spectrum of MRK 421 obtained 
in 1994 with the SW spectrometer. The best fit 
power law model is overplotted (solid line). The 
lower panel shows the residuals to the fit. 

Fig. 2. — x 2 contour confidence levels for normal- 
ization and energy index of the 1994 spectrum. 
The contours correspond to x 2 1lin +l 1 +4.61 and 
+5.99, respectively the la uncertainty for each 
single parameter and the 90% and 95% joint con- 
fidence levels. The best fit is marked by the filled 
symbol. 

Fig. 3.— April 25 - May 13 1995 Deep Survey 
lightcurve binned over one average EUVE orbit 
(~ 5544 s). The "flare", "decay" and "quiescence" 
intervals are marked. 

Fig. 4.— EUV spectrum of MRK 421 obtained in 
1995 with the SW spectrometer during a 0.3° off- 
axis observation. Hence note the increased cov- 
erage toward shorter wavelengths. The best fit 
power law model (dotted line) and the best fit 
power law + inverted gaussian model (solid line) 
are shown together with the residuals to the fits. 

Fig. 5. — Same as Fig. 2 for the 1995 spectrum 
and the power law model. 

Fig. 6. — Individual sections of the 1995 spectrum 
of MRK 421 (a) During the flare (MJD 49832- 
49835 or Apr. 25-28); (b) during April 29 - May 
6(MJD 49835-49843) for direct comparison with 
the analysis by Kartje et al. 1997; (c) during the 
remaining period, May 7-13 (MJD 49844-49851). 
The best fit power law model (dotted line) and 
the best fit power law + inverted gaussian model 
(solid line) are shown. The open circles in (b) 
correspond to the spectrum corrected with the on- 
axis effective area (not accurate for this off-axis 
observation) as presented in Kartje et al. 1997. 

Fig. 7. — Same as in Fig. 2 for the three sections 
of the 1995 observation and the power law model. 
They are displayed, from the top to the bottom, 
in temporal order: part 1 corresponds to the flare, 
part 2 to the decay and part 3 to the quiescent 
period. 

Fig. 8.— EUV spectrum of MRK 421 obtained 
in 1996 with the SW spectrometer during a 0.3° 
off-axis observation. The best fit power law model 



(dotted line) and the best fit power law + inverted 
gaussian model (solid line) arc shown. 

Fig. 9. — Same as in Fig. 2 for the 1996 spectrum 
and the power law model 

Fig. 10.— ASCA spectrum of MRK 421 obtained 
in 1995 during the flare. The solid lines represent 
the best-fit model when an absorbed power-law 
with a gradually changing spectral index is applied 
to the spectrum of all the instruments simultane- 
ously (upper panel ). The lower panel shows the 
residuals. 

Fig. 11. — EUV source flux vs. energy index (a) 
for the power law component of the power law plus 
gaussian model. 
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Table 1 

EUVE Observations of Mrk 421 



Year 


Instruments 


Start Date 


End Date 


rp b 

-L exp 


Count rate c 


Spectral Data 










(ks) 


(c s- 1 ) 




1994 d 


DS/S 


April 2 


April 12 


280 


0.18 


yes 


1995 


Scanner 


Feb. 4 


Feb. 7 


68 


0.45 


no 


1995 


DS/S 


April 25 


May 13 


355 


0.29 


yes 


1995 (flare) 


DS/S 


April 25 


April 28 


86 


0.37 


yes 


1995 (decay) 


DS/S 


April 29 


May 6 


154 


0.29 


yes 


1995 (quiescence) 


DS/S 


May 7 


May 13 


115 


0.22 


yes 


1996 


DS/S 


April 17 


April 30 


299 


0.30 


yes 


1996 


DS/S 


May 10 


May 11 


3.6 


0.30 


no 


1997 


Scanner 


Feb. 7 


Feb. 11 


108 


0.27 


no 


1998 f 


DS/S 


April 19 


May 1 









a "DS/S" indicates that the source was observed simultaneously in the 3 EUVE spectrometers and in the Deep 
Survey photometer. "Scanner" indicates that the source was observed simultaneously in the three EUVE photome- 
ters during a pointing within the Right Angle Program. 

b Total exposure time calculated eliminating (i) all SAA passages (ii) all satellite daytime data (iii) all times during 
which the detector was turned off (iv) all times affected by possible earth blockage. It does include corrections for 
telescope vignetting, deadtime and limited telemetry allocation (primbshing) . 

c Average count rate in the Lexan/B filter (« 60 — 180 A). The count rate in the scanner has been normalized 
to the one in the DS assuming an average spectral shape, to take into account the difference in the effective area 
between the two instruments. 

d This observation was presented in Fruscione et al. (1996 

e This part of the observation was presented in Kartje et al. 1997) 

f This observation is not included in this paper 
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Table 2 

Spectral fits to the EUV spectrum of MRK 421 assuming an absorbed power law model. 
Fixed parameters are N h = 1.45 x 10 20 cm~ 2 , N HeI /N HI =0.1 and N HeII /N HI =0.01 



Obs. 






Range 


a 


J 80 A 










(A) 




(phot cm~ 2 s _1 A -1 ) 




1994 






75-110 


0.85±°f 2 


6.44 ±0.33 


1.06 


1995 






70-110 


1.13 ±0.38 


6.70 ±0.22 


2.45 


Flare 


Apr. 


25-28 


70-110 


1.43 ±0.56 


10.28 ±0.51 


1.57 


Decay 


Apr. 


29-May 6 


70-110 


2.19 ±0.50 


7.24 ±0.33 


1.24 


Quiescence 


May 


7-13 


70-110 


-0.60±iil 


3.53 ±0.35 


1.47 


1996 






70-110 


2.12 ±0.31 


9.32 ±0.27 


1.61 
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Table 3 

Spectral fits to the EUV spectrum of MRK 421 assuming an absorbed power law model plus a reverted gaussian. Fixed 

PARAMETERS AS IN TABLE 2 



— 3 phot \ 



1994 -0.29 ± 0.83 6.72 ± 0.35 79.3±J;g| 0.20 ± 0.05 -10.67 ± 0.95 1.14 5.06 a 

1995 1.14 ±0.41 7.55 ± 0.22 73.5 ± 0.84 6.07 ± 1.47 -2.23 ±0.48 1.17 8.86 
Flare Apr. 25-28 0.71 ± 0.61 11.6 ± 0.52 74.87 ±0.45 3.04 ± 0.95 -5.86 ±0.14 0.57 14.11 
Decay Apr. 29-May 6 1.62 ±0.59 7.84 ± 0.34 72.75 ±0.28 1.86 ±0.51 -5.22 ± 1.30 1.17 9.33 
Quiescence May 7-13 0.63 ± 1.29 3.90 ± 0.35 71.43 ±0.04 0.17±^ -6.55 ±3.17 1.29 4.83 

1996 -0.15 ± 0.38 9.64 ±0.28 71.82±g;^ 0.66 ±0.12 -10.40 ±2.75 1.10 12.02 



a This is calculated in the 75 — 110 A range for the model in Table 2 



